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INTEGRATION OF THE RELATIVISTIC EQUATIONS OF 
MOTION OF AN ARTIFICIAL EARTH SATELLITE 

by 

Abolghassem Ghaffari 
ABSTRACT 

The Lindstedt perturbation method is applied to the 
motion of an artificial earth satellite which moves along 
a geodesic of the Schwarzschild metric of general rela- 
tivity. The purpose of this analysis is to determine the 
extent to which general relativistic effects are detectable 
in range measurements of earth orbiting spacecraft. 
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SUMMARY 

The Lindstedt perturbation method is applied to the 
motion of an artificial earth satellite which moves along 
a geodesic of the Schwarzschild metric of general rela- 
tivity theory. It is shown that both the amplitude and the 
frequency of the first approximate solution obtained is 
affected by the nonlinearity of the relativistic term ap- 
pearing in the equations of motion. The approximate 
periodic solution is compared with the solution for the 
motion of an artificial earth 24 hours synchronous satel- 
lite orbit (Application Technology Satellite 3) in a New- 
tonian force field. Assuming that the coordinate time and 
the initial conditions were the same in both systems, it is 
deduced that the maximum of the magnitude of the radial 
deviation in both systems is of order 1.6 cm after half 
orbit period and then too small to be detected. 


INTEGRATION OF THE RELATIVISTIC EQUATIONS OF 
MOTION OF AN ARTIFICIAL EARTH SATELLITE* 

1. Introduction. In the theory of general relativity, the external gravitational 
field of a spherically symmetric massive body M, whose center lies at r » 0, is 
represented by the static Schwarzschild metric 

ds 2 

An exact solution, for such a metric, of a set of discrete bodies is possible 
only when one of the bodies is of finite mass whereas the rest are of infinitesimal 
small mass [ 1 ] f 

If the massive body is taken to be the Earth, then an artificial earth satel- 
lite may be considered as an infinitesimal test particle whose gravitational field 
can be neglected. 

The complete solution for the relativistic effects of the combined fields of 
the Earth and the Sun on an artificial earth satellite is complex. The metric 
must involve both the principal masses, and the Earth is orbiting the Sun, the 
field is no longer static. It is well known [2] that in general the motion of a 

*Some results of this paper have been presented at the Sixth Semi-Annual Astrodynamic Conference 
held at Goddard Space Flight Center, Greenbelt, Maryland on November 7 and 8, 1967. 

^Figures in brackets indicate the literature references at the end of this paper. 
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dr 2 + r 2 (d<9 2 + sin 2 q d</ 2 ) I. 
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satellite can be derived with a great deal of precision by considering the effects 
of the Earth and the Sun separately. 


The relativistic effects of solar gravitation have been calculated by 
Corinaldesi and Papapetrou [3] and Papapetrou i4j. 


The use of artificial earth satellites seems more suitable than solar orbiting 
satellites for testing the effects of general relativity for the Sun has many 
planets and their mutual attractions are significant. 


For convenience we suppose the earth gravitational field is spherically sym- 
metric, and the Sun’s gravitational effect on an artificial earth satellite is negli- 
gible. Bearing in mind the above mentioned assumptions, one knows that the motion 


of an artificial earth satellite is given by the equations of the ordinary geodesics of 
the Schwarzschild metric [X, 6] which in simplified form is 


r 

d«jfc 


2 — L. + 2m •£— -i - y 

r 3 r 2 h 2 r 


or 


^du\ 2 0 3 2 2mc 2 

] = 2 mu ' 5 - u* + u - y 

h2 


( 2 ) 


where 


u -• 


•X - — (1-/3 2 ), 
h 2 


n- i-z 


- 1/2 


m 


GM 


h = r 


2 d0_ 
ds 
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r is the distance of the satellite from the center of the Earth, M is the mass of 
the Earth, v is the speed of satellite and e indicates the local speed of light, 
h and , are constants of integration, the angle / is the true anomaly and G 
is the gravitational constant. 

2. Integration of the Equations of Motion 

The classical equation (2) has been investigated and discussed completely in 
most star Hard texts on celestial mechanics and general relativity [5, 6, 7, 8 and 9). 
The rigorous integration of (2) leads to the Weirstrass elliptic function a (4, g 2 , g 3 ) 
which satisfies the differential equation 



4 o 3 - g 2 o - g 3 , 


where g 2 and g 3 are constants. However, in practice an approximate integration 
of (2) gives us the advance of the perigee. 

The solution of equation (2) with reference to the Sun as a central body, and a 
planet as the test particle gave the perihilion advances in 100 years of the planets 
Mercury (43?03), Venus (8 f . f 64), Earth (3 , . , 84) and Mars (lV35), predicted by 
Einstein’s law of gravitation, 

A complete discussion of equation (2) lias been given by J. Chazy [5] , 
and McVittie [6] discussed and solved equation (2) by a different method. 
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Differentiating Eq. (2) with respect to the true anomaly / and setting 
m c 2 /h 2 -v l/p, we get 

d 1 <> 2 

* ii — . * 3mu * . /m 

<1/ 2 P 

provided du/<W has only isolated zeros, i.e., eliminating the circular orbital 
solutions. 

In the Newtonian law of gravitation the right-iiand sides of (2) and (3) are 
quadratic and linear in u respectively. The presence of corrective terms 2 mu 3 
in (2) and 3 mu 2 in (3) is due to the Einstein’s law of gravitation. 

Equation (3) was first derived by A. Eddington [7] and using a method 
of successive approximation he obtained, in ignoring the small term 
3 mu 2 , as a first approximation a Keplorian orbit, then substituting the first 
approximation in the small term 3 mu 2 he arrived at the second approximation 
with a secular term which presents the resonance case. P. Bergmann [9] 
also considered equation (3) and applying Fourier Series procedure he ob- 
tained the perihilion advances of the planets without giving the explicit expres- 
sions of planets orbits. 

The purposes of this paper are: 

1. Integration of equation (3) by the classical Lindstedt method, and obtain 
a first approximate periodic solution of (3) starting from the perigee. 
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2. Comparison of the approximate periodic solution with iho solution for 
the motion of an artificial earth satellite in a Newtonian force field, 

3, In order to have a numerical estimate of the comparison the central body 
is taken to bo the Earth, and the test particle being identified with a 
long life earth synchronous satellite such as the Application Technology 
satellite 3. 

The advantages of equation (3) over equation (2) are twofold; 

a. Though equation (3) is of the second order, nevertheless it is nonlinear 
in u only, whereas equation (2) is nonlinear in both u and du/dv and its 
integration leads to the Weierstrass elliptic functions. 

b. We are primarily concerned with the periodic solution and its behavior 
of periodicity. The left-hand side of equation (3) suggests clearly the 
application of classical and modern approximation methods for derivation 
of approximations to periodic solutions. 

3. Appli cation of the Lindstedt Method. Let us assume that we wish to find an 
approximate periodic solution of Eq. (3) starting from the perigee i.e., satisfying 
the initial conditions: 


u(0) 


= =o • 

p \dtp/ 




In order to find such a solution and investigate the perturbation in the basic 


frequency arising from the presence of the relativistic term on the right-hand side 
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of equation (3), and also eliminate the secular term, we apply [10 3 the Undated! 
method to equation (3), 

The Lindstedt method consists primarily of changing the independent variable 
<p to another independent variable a such that the determination of the available: 
unknown coefficients enables us to eliminate gradually the secular terms in the 
subsequent approximations. As the frequency is altered it is of advantage fill 
to replace the independent variable 0 by a new independent variable a through 
the relation: 

^ = a (1 + Cj e + c 2 e 2 . . . ), (5) 


where o lf c 2 are unknown coefficients and e is an arbitrary small positive 

parameter. The smallness of the gravitational radius m enables us to 
suppose e = 3m, and then the equation (3) becomes 


or 


d 2 u 

da2 


d 2 u 

da 2 


(l + Cj e + c 2 e 2 , , ,) 2 ^u-i-eu 2 ^ = 0 

u -i + e ^ ~ + 2cj u - u 2 j + 0(e 2 ) = 0 


(6) 


Now let us write the solution of (6) in a power series with respect to the small 
parameter e : 


00 

u(e, a) = / e n u n (a), 

mwmmI 

n » 0 


( 7 ) 






and limit ourselves to the first order approximation such that 

u - u 0 + t Uj i 0 (' 2 ) . (8) 

Using the equations (G) and (7), one finds that the leading term u 0 is a solu- 
tion of the unperturbed equation 


d2ll 0 

da 2 


+ U Q 



( 9 ) 


\, 


and Uj satisfies 

d 2 u, o , 2c, 2o c 2 

i + u. = f 1 ° ~ 1 + — - (1 - c. p) cos a + cos 2a. (10) 

da 2 2p 2 p P 2 2p 2 

The unknown coefficient c t is to be determined such that no secular term appears 
in the solution of (10). Hence we chose c so that 



and therefore equation (10) assumes the form 


( 11 ) 


d 2 u. i 

1+ u. ~ (e 2 - 2 + e 2 cos 2a) 

da 2 2p 2 


( 12 ) 


which has the solution 


u 


i 2 - 2 


2p 


( 1 -cos 2 a) 


(13) 
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satisfying the initial conditions 


u j (0) 



0. 


( 14 ) 


Therefore the first order approximate periodic solution of (3), satisfying the 
initial condition (4), is 


U “ U 0 4 tUj + 0O- 2 ) 

where 


1 + ecos a 3 m(e 2 - 2 ) 

4 (1 „ cos 2a)+ 0( e 2 ) , (i 5 ) 

2 jp 


a 


1 + c 


. r r-7— + °( £2 ) = f 1 ~ t ) * f 0(t ' J) • < 16 > 


We can deduce that the nonlinearity shown by the relativistic term in (3) affects 
not only the amplitude of the solution but also the frequency. 

Formulas (15) and (16) show that the change in the frequency depends upon 
the amplitude a and the eccentricity e of the Keplerian orbit and also of the 
parameter e = 3 m, a property which belongs to the periodic solutions of all non 
linear autonomous differential equations. 

The period of the approximate periodic solution of (9) is equal to 2w , that 
is to say, the orbits are closed. The period of the exact solution of (3) differs 
from 2 77 by a small amount 8 , which is the difference between the angles oi two 
succeeding perigees, given by 


8 


"* 


2rr 

j 3m 
"p” 


- 2u 


2 'IT 



6m Tt 

P 


6nm 

a(l - o 2 ) 


rad. per revol. 


(17) 


Therefore the precession of the perigee of satellite orbit obtained by this method 
amounts to 6m Vp radians per revolution, which is in good agreement with the 
precession predicted by Einstein’s theory of general relativity as well as with 
gravitational theories of Whitehead [12] and Birkhoff [13], 


4. Comparison with Newtonian Force Field and Numerical Analysis 

Now let a second artificial earth satellite with the same characteristics 
features as the first one, move in a planar elliptic orbit according to Newtonian 
law of motion. If r N and <P are the classical polar coordinates of the second 
satellite in orbital plane, and u N - l/r N , the classical Newtonian equation of 
motion is 


u, 


1 + e cos (p 
P 


(18) 


where the true anomaly is measured from the perigee and 

p s a (1 - e 2 ) 

is the semi-latus rectum, a and e are the semi-major axis and the eccentricity 
of the Newtonian orbit respectively. We assume both satellites start from the 
perigee, i.e., they satisfy both the initial conditions: 
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f 


( 19 ) 


( r N ) 0 

< 1 , n)o 


( r R)o 

( u r)o 


> 


( U n) 0 


K) 


R' 0 




where u R and u N stand for relativistic and Newtonian solutions respectively. 

We are now in a position to compare the approximate periodic relativistic 
solution (15) with the classical Newtonian solution (18) by forming the difference 



The difference Au is a function of the Keplerian orbital elements and the gravi- 
tational radius m of the central body. As the primary body is taken to be the 
Earth, then m = 0.443 cm and the parameters a and e are constant for one re- 
volution. Therefore, the difference Au (<fi) is a function of the true anomaly <i> 
only for one revolution. 

In order to have a numerical estimate of Au (<£), formula (20) was applied to 
the case of the Earth as the central body and the Application Technology Satellite 
3 with the elements: 


a = 6.6109161 x R 
e = 0,0001703684 
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as a test particle, R ~ (>.371 x lo 8 cm is the mean radius of the Earth. The 
results are plotted and shown in Table 1 and Figure 1 in the range 0 ~ / - 180° 
and the function Au( 7 ), being a continuous function of •/ , begins to increase 
again after 180°. 

Now setting 


r R r r N «* a ^ r N 



( 21 ) 


where a- denotes the relativistic correction in radial deviation in both systems, 
one finds that 


Au = u N - u R 




( 22 ) 


Table 1 and Figure 1 show that, in case of Application Technology Satellite 3, 
the relativistic correction to equation (3) is too small to be detected and the 
maximum value of A u ~ 1.5 x 1CT 19 is obtained for <p = 89?59 , 59 , . , 148792708, 
correct to the nine decimal places given in the seconds. Therefore the maximum 
of the radial deviation in both systems* is of the order 

a ^ r 2 (1.5 x 1CT 19 ) a (4 x 10 9 cm) 2 (1.5 x 10” 19 ) ^ 1.6 cm , (23) 

and then too small to be detected. However, it is important to separate the 
relativistic effects from those due to the other causes, such as the oblateness 
of the Earth, magnetic fields, and the influence of high-altitude winds, etc. 
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Table 1. Variation of Au in Term? of the True Anomaly 4 > 
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Application of the averaging method of Bogoliubov-Mitropolsky 114], 
which leads to a different approximate solution, and its comparison with the 
solution (15) should follow later. 


5. Aclmowledgment 

I would like to thank Mr. Paul E. Schmid, Jr. for reading the manuscript and 
for many useful suggestions. 



6 . 


References 


1. MeVittie, G. C., Relativistic effects on Space Tracking Data. First Semi- 
Annual Status Report, NASA NGR 14-005-088 (1907), to appear in the 
Astronomical Journal. 

2. Sitter, W. de., Monthly Notices, Roy. Astron. Soc. 77, 155 (1916), 

3. Corinaldesi, E. and Fapapetrou, A. Proc. Roy. Soc. A. 209 (1951) 259. 

4. Papapetrou, A., Proc. Roy. Soc. A-209 (1951) 248. 

5. Chazy, J. La th^orie de la relativite et la m£canique Celeste, Tames I 
et II, Gauthier-Villars, Paris (1928) et (1930), 

6. MeVittie G., General relativity and Cosmology, 2nd Ed. Illinois University 
Press, Urbana (1965). 

Ml 

7. Eddington, A., The mathematical theory of relativity, Cambridge University 
Press (1924). 

8. Synge, J. L., Relativity: The general theory, 2nd Printing (1964), North- 
Holland Publishing Company, Amsterdam. 

9. Bergamn, P., Introduction to the Theory of Relativity, Prentice-Hall, Inc. 
(1942). 


15 


10. Ghaffari. A., On Home Applications of Approximate Methods in Relativistic 
Celestial Mechanics. To appear in the Proceedings of the Fifth International 
Conference uu Nonlinear Oscillations, Kiev USSR, August 25- .September 5, 1969. 

11. Bellman, R. K., Perturbation Techniques in Mathematics, Physics, and 
Engineering, Ilolt, Rinehart and Winston, Inc. (3904) . 

12. Whitehead, A. N., The principle of relativity, with applications to physical 
Science, Cambridge University Press (1922). 

13. Birldioff , G, D. , Relativity and Modem Physics, with the Corporation of H, K. 
Langer, Harvard University Press, Cambridge (1923). 

14. Bogoliubov-Mitropolsky, Asymptotic methods in the Theory of Nonlinear 
Oscillations. Translated from Russian (revised ed. 1958) by Gordon and 


Breach (1962). 


